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Abstract 



We point out that the suppression in the ratio i?dAu recently observed by the Brahms 
collaboration in forward scattering is stronger than usually appreciated. This is related 
to the fact that at forward rapidities Brahms measures negatively charged hadrons and 
that -RdAu is defined from the ratio of dA and pp scattering cross sections. We also 
investigate the influence of standard shadowing on i?<jAu) anc l the- typical values of partonic 
momentum fractions relevant in forward scattering. We find that iau > 0.02 dominate in 
the cross section. 



I Introduction 



The Brahms collaboration has recently presented measurements of the ratio -RdAu of hadron 
production cross sections in dAu and pp collisions pp. With increasing hadron rapidity, the 
data indicate a growing suppression of the ratio. Theoretical analyses of the data have focused 
on nuclear effects related to the gold nucleus, emphasizing variously the role of small x physics 
in the nuclear production [2], as well as of nuclear-enhanced power corrections in the dA cross 
section [3]. Other studies [HE] have addressed the question of whether leading-twist shadowing 
could be responsible for the observed suppression. 

In the present paper we hope to add valuable information that will have an impact on the 
present and future attempts to interpret the Brahms data, and on plans for further measure- 
ments. We base our analysis on next-to-leading order (NLO) leading-twist perturbative-QCD 
(pQCD) calculations of inclusive hadron production. Such calculations have enjoyed consid- 
erable success in comparisons with data from pp collisions at RHIC at -^i = 200 GeV. They 
yield good agreement with the Phenix data |Hj for pp — > tt°X at central rapidities, and with 
data from Star j2| on pp — > tt°X at forward rapidities r] = 3.8 and rj = 3.3, the latter being 
essentially in the kinematic range explored by Brahms in their most forward measurements. It 
is appropriate to point out that there are sizable uncertainties in the NLO calculation, related 
to the choice of fragmentation functions and scales. However, we are confident that for the 
kinematics relevant for Brahms NLO pQCD does explain at least 50% of the cross section, 
as a conservative estimate, and hence is a viable tool for obtaining deeper insights into the 
production mechanism, at least for pp collisions. 

The main point of our analysis can be stated very simply: the nuclear effect reported by 
Brahms actually appears to be substantially larger than appreciated in the studies 
This is related to the fact that in the very forward region, at rapidities rj = 2.2 and rj = 3.2, 
Brahms only measured negatively chargedhadrons (h~) and not the charge average (h + + h~) /2 
as at the more central rapidities. In the ratio RdA the denominator refers to pp collisions, and 
negatively charged hadrons are expected to be produced more rarely in pp than in dA collisions, 
for which from isospin considerations it follows that at least 7r + and 7r~ should be produced in 
equal numbers. This immediately implies that the ratio of dAu — > h~X and pp — > h~X cross 
sections is "intrinsically enhanced", by actually a factor of about 1.5, as we will show. The fact 
that Brahms nonetheless reports a suppression of the ratio is therefore truly remarkable and 
awaits further investigation. 

We also revisit in our analysis the following questions: 
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(i) what are the most relevant parton momentum fractions x for hadron production at BRAHMS, 

in particular at very forward rapidity 77 = 3.2 where the suppression of R&x is largest? 
To what extent are truly small x, say, x < 10~ 3 probed? 

(ii) how relevant is leading-twist nuclear shadowing for the explanation of the Brahms data? 

These questions have already been addressed in some detail in pj. Our analysis extends that 
study by providing results within a full NLO calculation. This will generally lead to more 
reliable results. In addition, the enhancement effect mentioned above has of course also direct 
implications for estimates for R^a obtained when using leading-twist nuclear shadowing. Our 
calculations therefore provide an improved estimate as compared to the results of [5 , where 
the enhancement was not taken into account. We furthermore explore more thoroughly the 
possible effects of leading-twist nuclear shadowing by making more extreme assumptions on 
the structure of screening at intermediate x, tolerated because of possible uncertainties in the 
connection between diffractive HERA data and gluon shadowing. 

In sec. iniwe discuss the ranges of partonic momentum fractions mainly probed by the for- 
ward Brahms data. Section ITTTl addresses the leading-twist nuclear shadowing and its relevance 
in forward dA scattering. With the findings of sections |H] and IIII1 we are in the position to 
discuss RdA in the forward region in more detail. This is done in Section HV| where we empha- 
size our main point related to the normalization of RdA by the pp — > h~X cross section. We 
summarize and conclude with section IVl 



II Kinematics and x ranges probed in forward scattering 

We consider the reaction H1H2 — > hX, where Hi, H% are initial hadrons and h is a hadron in the 
final state produced at high transverse momentum p?. Since large pt ensures large momentum 
transfer, the cross section for the process may be written in a factorized form, 

da = do ° 2 dxi / dz fa 1 ^^) fb 2 ( x ^^) D c( z ^) 



xdal b (xiP Hl ,x 2 PH 2 ,Ph/z,n) , (1) 

where the sum is over all contributing partonic channels a + b — > c+. . ., with do c ah the associated 
short- distance cross section which may be evaluated in QCD perturbation theory: 

dK b = daT + + • • • " ( 2 ) 

7T 
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The leading-order (LO) contributions da c Jp are of order a 2 s \ the next-to-leading order (NLO) 
corrections are known jH] and will be included in our analysis. 

In Eq. (JTJ), f^i.x.n) denotes the distribution function at scale \i for a parton of type i in 
hadron H, carrying the fraction x of the hadron's light-cone momentum. Likewise, D^(z,fi) 
describes the fragmentation of produced parton c into the observed hadron h, the latter taking 
momentum fraction z of the parton momentum. The scale \x in Eq. stands generically for 
the involved renormalization and factorization scales. \x should be of the order of the hard scale 
in the process; in the following we choose fj, — Pt- The dependence on /i is actually quite large 
even at NLO [B] ; however, in this work we are mainly interested in ratios of cross sections for 
which the \x dependence is fairly insignificant. 

The lower limits of the integrations over momentum fractions in Eq. (pQ) may be derived in 
terms of xt = 2px/\/s and the pseudorapidity 77 of the produced hadron. They are given by 
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From these equations it follows that at central rapidities rj « the momentum fractions x\ and 
X2 can become as small as roughly Pr/y/s. In forward scattering, that is, at (large) positive rj, 
the collisions become very asymmetric. In particular, X2 may become fairly small, whereas x\ 
tends to be large. For forward kinematics at Brahms one has, typically, Pt ~ 1.5 GeV and 
rj = 3.2. This implies that X2 may become as small as ~ 3.5 x 10 -4 . However, in practice it 
turns out that such small X2 hardly ever contribute to the cross section: if X2 is so small, the 
hadron with transverse momentum can only be produced if both X\ and z are unity, where 
however the parton distributions f^fai,^) and the fragmentation functions D^(z,fi) vanish. 
This is an immediate consequence of kinematics, as demonstrated by Eqs. 0- One can show 
that if the parton density f^ 1 (xi,fi) behaves at large x± as (1 — xi) af and D^(z, (i) as (1 — z) aD 
(with some powers af,a^ 3> 1), the ^-integrand in Eq. (JTJ vanishes in the vicinity of x™ n as 
(x2 — x 1 £ m ) a f +aD+1 . Therefore, contributions from very small X2 are highly suppressed. 

The question, then, remains of how small X2 really is on average for forward kinematics 
at RHIC. This is of course relevant for judging various explanations for the suppression of 
i?dA seen by Brahms, in particular those relating to saturation effects in the nucleus wave 
function [2] . Figure^shows the distribution of the cross section for pp — > tc°X at y/s = 200 GeV, 
Pt = 1-5 GeV, rj = 3.2, in bins of log 10 (x2). The overall normalization is unimportant of 
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Figure 1: Distribution in log 10 (a;2) of the NLO invariant cross section Ed 3 a/dp 3 at y/s = 200 GeV, 
= 1.5 GeV and rj = 3.2. 

course; for definiteness we note that the sum of all entries shown in the plot yields the full 
NLO invariant cross section Ed 3 a/dp 3 in pb/GeV 2 . For the calculation we have chosen the 
CTEQ6M jH] parton distribution functions and the fragmentation functions of Ref. ^U]. One 
can see that the distribution peaks at x 2 > 0.01. There are several ways to estimate an average 
(x 2 ) of the distribution. For example, one may define (x 2 ) in the standard way from evaluating 
the integral in Eq. $IJ with an extra factor x 2 in the integrand, divided by the integral itself: 

f x min dX2 X2 iff 2 (£2, /i) • • • 



fa) 



(4) 



J xT n dx 2 f^ 2 (x 2 ,fi) . . . 

where the ellipses denote the remaining factors in Eq. ((TJ). Alternatively, one may simply 
determine (x 2 ) as the median of the distribution, demanding that the area under the distribution 
in Fig.^to the left of (x 2 ) equals that to the right. Either way, one finds an average (x 2 ) > 0.01, 
typically 0.03 — 0.05 at this Pt and 77. 

The precise shape of the distribution and the value of (x 2 ) depend somewhat on the parton 
distributions (and, less so, on the fragmentation functions) chosen. We remind the reader that 
the distribution shown in Fig. ^is at p? — 1.5 GeV and that we have chosen the factorization 
and renormalization scales to be \i = pr- This means that we are using a fairly low scale in the 
parton densities. At this scale, the CTEQ6 densities, in particular the gluon, are still relatively 
flat towards small x. In order to estimate to what extent this influences the distribution, we 
have calculated it for the GRV jTTj parton distributions, which are steeper at this scale. The 
corresponding histogram in log 10 (x2) is shown in Fig. 121 One can see that as expected it peaks 
somewhat more to the left; nevertheless there is not much quantitative change in the average 
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Figure 2: Same as Fig. [T] but for the GRV [llj parton distributions. 



x 2 . The full invariant cross section is about 15% smaller than for the CTEQ6 set. We have 
mentioned in the introduction that there are data from Star for the cross section for pp — > 7r°X 
in roughly this kinematic range [Zj which are in very good agreement with the NLO calculation 
used in Fig.^ This supports the view that the distributions shown in Figs. ^ and 121 are realistic. 



Figure El shows the log 10 (x2) distribution at px = 5 GeV. At this pt, one is closer to the 
boundary of phase space given by the condition xrcosh^) = 1, where all momentum fractions 
x\, #2) z are forced to 1. The distribution in X2 is therefore more "squeezed" and shifted to the 
right. The effect is countered to some extent by evolution since at scale /i = 5 GeV the parton 
distributions are steeper than at /x = 1.5 GeV. 

Finally, in Figure |U we present results for the averages of x\, x 2 , and z, as functions of 
pion transverse momentum and rapidity at ^/i = 200 GeV. Here we have defined (x 2 ) as in 
Eq. (jlj), with analogous definitions for (x\) and (z). The upper row shows results for fixed 
Pt in forward scattering. Besides rj = 3.2 as relevant for Brahms, we have also extended the 
results to rj = 4.2 which may be useful for future experimental studies. It becomes evident 
that (xx) and (z) are very large in forward scattering, as expected, (z) is particularly large 
because, on account of Eq. the single fragmentation function has to compete with two 
parton densities, each function being strongly suppressed at large momentum fraction. As we 
have already seen in the histograms, Figs. (^2) is much smaller. As px increases and the 
boundary of phase space is approached, all momentum fractions become larger and eventually 
converge to unity. We also note an unexpected upturn of (x 2 ) toward smaller p T . We have 
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Figure 3: Same as Fig. [U but for p T = 5 GeV. 



not been able to identify this effect as resulting from any straightforward origin. The precise 
small-x behavior of the parton distributions at the rather low scales involved here plays a role 
(however, the effect also occurs for the steeper GRV distributions). The structure of the cross 
section formula in Eq. ((J) itself is also partly responsible. In the lower part of Fig. 0] we show 
the averages as functions of rapidity for two fixed values of pt- At rj = one obviously starts 
from (xi) = (#2); with increasing rj the two momentum fractions become very different. Toward 
77 = cosh - (l/a>r) they again both tend to unity; for (x 2 ) this happens rather late. 



Ill Influence of leading-twist nuclear shadowing 

In dA collisions, nuclear effects will alter the distribution in x 2 a s well as the full cross section. 
Shadowing effects at small-a^, x 2 < 0.05 — 0.1, will lower the cross section and will lead to yet 
higher average x 2 being probed. There will also be enhancements at larger X2, 0.05 — 0.1 < x 2 < 
0.2, associated with yet another coherent nuclear effect, antishadowing. This will be followed 
by the suppression related to the EMC effect for 0.2 < x 2 < 0.8, and by the subsequent 
enhancement explained by the Fermi motion for x 2 > 0.8. 

Since in the Brahms kinematics the average x 2 is near 0.01, the principal nuclear effect is 
shadowing. We investigate its role in the interpretation of the Brahms data by considering 
leading-twist shadowing, using the parameterization of nuclear parton distribution functions 
(nPDFs) derived in ^21EIE]- Unlike most other sets of nPDFs [HE], these functions have 
a rather rapid and strong onset of shadowing toward small-x, so they may serve as a good tool 
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Figure 4: Averages of X\, x 2 , and z in pp — ► n°X at RHIC, defined as in Eq. (J4j), as functions 
of pion transverse momentum at forward rapidities (upper row), and of pion rapidity for fixed pr 
(lower row). 



for studying the "maximally possible" effects of leading-twist nuclear shadowing in forward dA 
scattering. The recent study jH] has also employed the nPDFs of fl2\ If"""} II"""] , albeit only in the 
framework of a lowest order (LO) calculation. 

We note that we will neglect any nuclear effects in the deuteron, for which we just use 
d = (p + n)/2. As follows from our analysis of the average x%, the deuteron parton distributions 
are sampled at values of x% in the domain of the EMC effect. Therefore, the approximation 
d — (p + n)/2 is valid to a few per cent accuracy, as can be estimated as follows. The CTEQ fits 
use the neutron structure function extracted from the deuteron data within the framework of 
the convolution approximation (Fermi motion). The convolution model gives for the structure 
function ratio R = 2F£/(F%+F£) the values of 0.99 for x ~ 0.5 and of 1.03-1.05 for x ~ 0.7. As 
a result, our d — (p + n)/2 approximation overestimates the true deuteron parton distributions 
by about ~ 1% at x ~ 0.5 and underestimates them by a few percent at x ~ 0.7. Since large xi 
are important in our calculations, proper account of this would make the effect we will discuss 
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in the next section even slightly bigger. Note that for heavier nuclei the convolution model 
contradicts the EMC effect. However here we are using it to "restore" the deuteron structure 
function within the procedure used to extract the neutron structure function; see fH] for an 
extensive discussion of nuclear effects in the deuteron parton densities. 

Let us now briefly describe the approach of ^21 El El for deriving nPDFs. Leading- 
twist nuclear shadowing is obtained using Gribov's theorem relating nuclear shadowing 
to diffraction, Collins's QCD factorization theorem for hard diffraction in DIS ^B], and the 
QCD analysis of hard diffraction measured at HERA in terms of diffractive parton distribution 
functions of the proton. Operationally, the nPDFs are first derived at the initial scale Qo=2 
GeV and for < x < 1. Standard (NLO) DGLAP evolution is then used to obtain the 
nPDFs for Q 2 > Qg. 

Analyses of DIS by both HI and ZEUS [20] demonstrate that diffraction constitutes 
~ 10% of the total cross section in the quark-dominated channel and as much as ~ 30% in the 
gluon channel. As a result, it is found in ^21 CSl E] that the effect of nuclear shadowing in 
nPDFs is large and, even more strikingly, much stronger in the gluon nPDF at small x than in 
the quark nPDFs. 

Conservation of the baryon number and the momentum sum rule then require that the 
depletion of nPDFs at small values of x, x < 0.01, be accompanied by a certain enhancement 
at 0.05 < x < 0.2. The transition from shadowing to enhancement, and the enhancement itself, 
are not described by the Gribov theorem and, hence, can be only modeled. Using the available 
fixed-target nuclear DIS data j2I] as a guide, the "standard" scenario of O CSl HI] assumes 
that the transition from nuclear shadowing to the enhancement takes place at x = 0.1 for quark 
nPDFs, and at x — 0.03 for the gluon nPDF. In the following, we will refer to this set of nPDFs 
as "shadowing 1". 

Figure shows the distribution of the NLO cross section for dAu — > 7r°X in log 10 (x2), using 
shadowing 1. For comparison, we also display the previous result for pp — > ir°X of Fig. E] (solid 
line). A clear shift in the distribution to larger x 2 is visible. At small X2, there are significant 
shadowing effects, while at large x<i there is a small contribution from antishadowing. It is 
evident from comparison of the areas underneath the distributions that the net effect on the 
dAu cross section will be a decrease. However, one can anticipate that the decrease will be 
rather moderate: while nuclear shadowing does deplete the dA cross section compared to the 
pp cross section, the probed values of x<i are clearly not small enough to deliver a significant 
nuclear shadowing effect. 
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Figure 5: Same as Fig. [U but also showing the result for dAu scattering using the shadowing 
of Refs. [T21 H3 HIT and a more extreme ansatz for shadowing (see text). The "spikes" in the 
distributions at log 10 ( a; 2) ~ —0.8 are artifacts of the implementation of antishadowing in the nPDFs 
of [121 E31 El- For better comparison we have displayed the result of Fig. [T] by a solid line. 

The kinematics for forward scattering at Brahms mostly corresponds to values of x^ in the 
transition region between shadowing and antishadowing, where the predictions for nPDFs are 
rather uncertain. Therefore, in addition to the standard scenario ("shadowing 1") of nuclear 
shadowing, we have also explored an option for which nuclear shadowing in the gluon channel is 
increased by extending it up to x = 0.1, similarly to the shadowing in the quark densities. We 
refer to the resulting set of nPDFs as "shadowing 2". The corresponding log 10 (x2) distribution 
of the cross section for dAu — > n°X is also displayed in fig. 03 Compared to shadowing 1, there is 
only a small modification of the distribution, which will lead to a very slight further suppression 
of the dA cross section. It is worth emphasizing that one can hardly increase the amount of 
gluon shadowing at x ~ 10~ 3 since here there are constraints from J/ip data [22] • The model 
of ^21E3E] gi yes a reasonable description of the observed suppression by a factor ~ 2, which 
would be spoilt by a much stronger gluon shadowing. At the same time, as soon as the amount 
of gluon shadowing at x ~ 10~ 3 is fixed, the gradual decrease of shadowing with increasing 
x is automatic as a consequence of the decrease of the coherence length oc Xjm^x and of a 
smaller probability for diffraction. Therefore, we conclude that the standard effect of leading 
twist nuclear shadowing will at best be able to explain only a small fraction of the dramatic 
suppression of the spectra of charged hadrons at forward rapidities observed by Brahms. This 
statement is in line with the LO result in the revised version of 

Figure El shows the corresponding results for p T = 5 GeV. Here, larger x 2 are probed, and 
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Figure 6: Same as Fig. 03 but at p T = 5 GeV. 

only slight antishadowing effects appear. 

So far we have only considered tt° production as an example. This is however not really 
appropriate for a comparison with the Brahms data which, at forward rapidities, are for 
negatively charged hadrons h~ . As we will now show, for the case of h~, even in presence of the 
shadowing effects just discussed, the leading-twist NLO calculation predicts an enhancement, 
rather than a suppression, of RdA- 

IV Isospin considerations for the ratio of dA and pp cross 
sections 

We now consider the ratio R^a of single-inclusive hadron cross sections in dA and pp scattering. 
The Brahms experiment has presented data pQ for R^a as a function of hadron transverse 
momentum p?, in four different bins of hadron pseudorapidity r), with central values r] = 
0, 1, 2.2, 3.2. Brahms sees a significant suppression of the ratio with increasing r\. 

While Brahms measures inclusive charged hadrons, (h + + h~) /2, at central rapidities (rj = 
and 1), their R^a data at forward rapidities refer only to negatively charged hadrons h~ . This 
has profound consequences. To see this, let us assume for the moment that pions dominate 
the spectrum of observed high-p^ hadrons. Negatively charged pions are produced more rarely 
than positively charged ones in pp collisions, due to the up-quark dominance in the proton. An 
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Figure 7: Ratio of pp — > n + X and pp — > 7r~X cross sections at -y/s = 45 GeV and scattering 
angle 50° (corresponding to pseudorapidity rj = 0.76), as a function of x T = 2p T /y/s. The data 
points are from measurements at the ISR |23j. The curve shows the result of the NLO calculation, 
using the fragmentation functions of |24j. 

example for this is shown in Fig.0 where we display data for the tt + /tt~ ratio from the ISR [23] 
at \/s = 45 GeV. We also show the result of the NLO calculation, using the fragmentation 
functions of Ref . [21] , which provides separate sets for negatively and positively charged pions. 
[We note that there are also tt + /tt~ data at ^/s = 62 GeV [213 which lie lower and are in less 
impressive agreement with NLO pQCD.] As Fig.[7|shows, there is clear excess of positive pions 
over negative ones. In contrast, isospin considerations imply that tt + and tt~ are produced 
practically equally in dAu collisions. Therefore, one expects R^a for negatively charged pions 
to be intrinsically enhanced, if it is normalized by the pp cross section and not, for example, 
by the dp one. 

To go into a little more quantitative detail, we recall that at forward rapidities the partonic 
collisions are very asymmetric. Large contributions to the scattering come from situations in 
which a large-xi valence quark in the "projectile" (i.e., in the deuteron, or in one of the protons) 
hits a small-a^ gluon in the "target" (i.e., in the gold nucleus or in the other proton). The 
underlying (LO) subprocess is then the quark-gluon Compton process qg — ► qg. For negatively 
charged pions one then expects that down quarks play a particularly important role in the 
Compton process, since both the "projectile" and the produced ir~ have a d valence quark. To 
a good approximation (see the previous section), the deuteron's d valence density is given by 

^deuteron _ (dF V(d + d™ al ) = - {dF val + U^ al ) , (5) 

where we have used isospin invariance to relate the valence-c? distribution in the neutron to the 
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Figure 8: Ratio d^ teron {x) / d p val {x) as a function of x at the scale ji = 2 GeV, as given by the 
CTEQ6M set. The deuteron's c/-valence distribution has been estimated according to Eq. (JHJ). 

valence-w in the proton. Due to the up-quark excess in the proton, the distribution in Eq. © 
becomes much larger than the proton's <i-valence distribution at high x, as shown in Fig. |SJ 
resulting in an enhancement in -RdA- Of course, other scattering channels will contribute as well 
and dilute this valence effect. In addition, Brahms does not measure only pions, but inclusive 
charged hadrons, = ir^ + K ± + P + . . .. Nevertheless, when changing to charged-hadron 
fragmentation functions as given by [21], we find that the difference in deuteron and proton 
valence densities continues to play an important role in the forward production of negatively 
charged hadrons. This is demonstrated by Fig. El The solid lines show the ratio R^a at rj = 
and 1 for summed charged hadrons (h + + h~) /2, and at r\ — 2.2 and 3.2 for negatively charged 
hadrons, exactly corresponding to the Brahms conditions. We have used the "shadowing 1" 
set described in the previous section. As shown in sec. [H] when going from r\ = to r\ = 1, 
the average x<i probed slightly decreases, and shadowing effects start to become visible at the 
smaller p?. At rj = 2.2 and r] = 3.2 the x<i become yet smaller, but since now negatively charged 
hadrons are measured, the valence effect discussed above outweighs any stronger shadowing, and 
in fact the ratio RdA strongly increases with pt because larger and larger x\ become relevant. 
For comparison we also show in Fig. El the results at r\ = 2.2 and rj = 3.2 for summed charged 
hadrons (dashed lines). For these, the effective valence densities in the proton and deuteron 
are the same, and the enhancement seen for negative hadrons disappears. Shadowing effects of 
up to 15% are visible then, as expected from Fig. 03 and as also found in Ref. [5] where only 
summed charged hadrons were considered. It is remarkable that for the highest px ~ 3 GeV at 
7] = 3.2 the curve for h~ is enhanced by about a factor 1.5 with respect to the one for summed 
charges. 
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Figure 9: Ratio R^a of cross sections for dAu — > hX and pp — > ftJT as a function of transverse 
momentum at various rapidities relevant to the Brahms experiment. As in experiment, we have 
considered production of summed charged hadrons, h ch - = [h + + h~)/2 for r\ — 0, 1 and negatively 
charged hadrons h~ for rj = 2.2 and 3.2. For comparison, the dashed lines show the result for 
summed charged hadrons at r] = 2.2 and 3.2. We have used the "shadowing 1" nPDFs for the gold 
nucleus. The fragmentation functions are from [24J; we have found that for the case of summed 
charged hadrons using the set of [10J does not alter our results by more than a few per cent. 

Our analysis so far has been entirely based on a NLO pQCD leading-twist calculation and 
the use of fragmentation functions extracted from e + e~ annihilation data. We expect that even 
if nonperturbative phenomena are important in the kinematic regime explored by Brahms, the 
enhancement in R^a resulting due the different "projectiles" in the numerator and denominator 
will persist. Here, our reasoning is as follows: ISR data (see [26J and references therein) on 
tt^ production at pr ~ 1 GeV and large Feynman-x^p ~ 0.4 were found to be consistent with 
da pp ^ 7T+ /dxp oc u(xp), da pp ^ 7T jdxp oc d(xp), where u and d are typical densities for up and 
down valence quarks. Such spectra in xf are much harder than the ones pQCD would give, 
indicating the presence of a nonperturbative production mechanism. Nonetheless, the ratio of 
7r + to 7T~ cross sections is of a similar magnitude as the one in the perturbative regime shown 
in Fig. [71 Adding the fact that the yields of n + , n° in dA scattering should be equal because 
of isospin, irrespective of the production mechanism, we are led to the conclusion that R^a for 
negatively charged hadrons should be enhanced even if nonperturbative effects dominate. Note 
here that nonperturbative effects in the fragmentation region are known to be consistent with 
the Feynman scaling. Hence if nonperturbative effects are important they should lead to about 
the same ir + /ir~ ratio in pp scattering for the same xp at ISR and at RHIC. 

We finally note that a potential caveat to our main finding in Fig. Incomes from a further set 
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of Brahms data. While we have mentioned that Brahms measures all hadrons and not just 
pions, we have generally assumed that pions dominate the observed hadron spectrum. However, 
preliminary data from Brahms 27. show that the cross section for dAu — ► h + X becomes 
significantly larger than that for dAu — > h~X at pt > 1 GeV. At pt = 3 GeV, they observe 
about three times as many h + as h~ . As we have pointed out before, isospin excludes that such 
an excess could be due to pions: cr(dAu — > n + X) = a(dAu — > tc~X), up to corrections of a few 
per cent related to the fact that the gold nucleus is not isoscalar. Standard sets of fragmentation 
functions do not predict large contributions from kaons and protons, and in the NLO calculation 
one ends up with <r(dAu — > h + X) at most only 10% larger than <r(dAu — > h~X). It is hard to 
conceive that proton production could lead to a large enhancement of h + over h~ (that even 
increases with p?) but, barring any experimental systematic problem, this appears to be the 
conclusion at present. If the final Brahms data continue to show this large excess, it will 
be a challenge to understand it in terms of a nonperturbative effect. Such an effect could, 
perhaps, result from coalescence of quarks from the incoming nucleon with other partons, to 
form a baryon. For this to work, the quarks would need to experience large transverse "kicks" 
and would need to lose a significant fraction of their momentum. Such a possibility could be 
connected to the expectations of a very strong suppression of the forward nucleon spectrum in 
central nucleon- nucleus collisions [2E1 EH] • 



V Conclusions and outlook 



We have shown that there is an intrinsic enhancement in the ratio -RdAu f° r negatively charged 
hadrons, simply because of the different nature of the "projectile" (deuteron vs. proton) in the 
numerator and denominator of -RdAu- In the light of Fig. El the significant suppression seen by 
Brahms at r\ — 2.2 and rj = 3.2 is even more striking than usually appreciated. The effect we 
have found has not been taken into account in any previous theoretical study [3 El E] of the 
data, to our knowledge. We expect that future data for -RdAu for summed charged hadrons in 
this kinematic regime will show an even stronger suppression than observed for h~ , roughly by 
a factor 1.5. The same should happen if -R p a u , rather than -RdAu, were measured for h~. 

Because of the effect, it is entirely impossible to explain the suppression in -RdAu at forward 
rapidities by a conventional modification of the leading-twist parton densities in nuclei. How- 
ever, even if we disregard the effect, nuclear leading-twist shadowing plays a rather unimportant 
role, giving at most a suppression of 15%. The reason for this is that parton momentum frac- 
tions in the gold nucleus are not very small on average even for forward kinematics, as we have 
shown. In other words, a large nuclear contribution originates from a range of x where nuclear 
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effects are known to be small (or even antishadowed). This generally sets severe limitations on 
the ability of any initial-state small-x effects to explain the observed suppression. 

We have also mentioned that the fact that Brahms observes the cross section for dAu — > 
h + X to be significantly larger than that for dAu — > h~X at pr > 1 GeV, indicates the presence 
of sizable nonperturbative contributions possibly related to protons. It appears likely that 
mechanisms responsible for an enhancement of proton production then also play a role in the 
observed suppression of -RdAu- Nonperturbative production of pions, too, could play a role: 
coalescence effects involving spectator partons are likely to be strongly suppressed in p(d)A 
collisions as compared to the pp case [28 . This suppression is further enhanced when energies 
are large enough to resolve the small- X high gluon densities [2H|- Such non-leading twist effects 
should decrease with increase of the transverse momentum of the pion, which is consistent with 
the trend in the data on -RdAu- On the other hand, the observed excess of dAu — > h + X over 
dAu — > h~X actually increases with p?, which is quite challenging to understand. Note also 
that the parton energy losses that would be necessary to reproduce the Brahms effect appear to 
be rather large: about 10% energy loss would be needed if it occurred only in the initial state 
before the interaction. If one assumes that the rate of energy loss is the same in the initial and 
final states, a loss of about 3% is necessary. For the kinematics relevant here, the suppression 
is more sensitive to losses in the final state since the average (z) for fragmentation was found 
to be substantially closer to one than the averages (#1,2) in the parton densities (see Fig. HJ). 

To further investigate experimentally the origin of the suppression in -RdAu it would be useful 
to perform measurements of dihadron production. Within LO kinematics, the pseudorapidities 
of the two hadrons are related by 771+772 = ln(xi/x2). One may therefore single out contributions 
from small £2 by demanding that both hadrons be rather forward. Fig. E3 shows this for 
a sample calculation. We assume that a "trigger" hadron (jr°) is detected with transverse 
momentum px,i — 2-5 and forward rapidity 2.5 < 771 < 3.5. Let the second 7r° have 1.5 GeV < 
Pt,2 < Pt,i- Without any restriction on the rapidity 772 of the second hadron, one then obtains 
the higher log 10 (a;2)-distributions in Fig. El As expected, these look very much like the single- 
hadron distributions shown in Sections|TTland lllTl If now the second hadron is also in the forward 
region at 1.5 < 77 2 < 4, the lower distribution is obtained, which is entirely located at Xi < 0.01. 
We also show in Fig. El the corresponding results for dAu collisions, using our "shadowing 1". 
The shadowing effects are much more relevant for the double-forward distribution, as expected. 
We note that the distributions in Fig. El are normalized such that they sum to the cross section 
da/dpT,i in pb/GeV. It is evident that there is a significant decrease in rate for two forward 
hadrons. The results shown in Fig. El are only LO. The NLO corrections are expected to be 
sizable; they are available [30 . The two hadrons we have studied in Fig. El would be nearly 
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Figure 10: LO distributions in log 10 (x 2 ) of the cross section da/dpT,i for pp — > 7r°7r°X and 
dAu — > TT°n°X production at y/s = 200 GeV. The kinematic variables have been chosen as described 
in the text. Solid lines are for collisions and histograms are for dAu, using "shadowing 1". The 
higher-lying histograms are for the case of "arbitrary" (i.e., unconstrained) r/ 2 , the lower ones are 

for 1.5 < 772 < 4. 

back-to-back in azimuthal angle. Further insights into the dynamics may be gained by studying 
back-to-back azimuthal correlations [HU IH2] • We finally note that small-x 2 effects might also 
become more readily accessible in conceivable future pA collisions after increase of the proton 
energy to 250 GeV at a later stage of RHIC operations. 
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